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Abstract We present a multi-line study of the massive star-forming region IRAS 
22506+5944. A new 6.7 GHz methanol maser was detected. 12 CO, 13 CO, C I8 and HCO + J 
= 1-0 transition observations reveal a star formation complex consisting mainly of two cores. 
The dominant core has a mass of more than 200 M , while another one only about 35 M . 
Both cores are obviously at different evolutionary stages. A 12 CO energetic bipolar outflow 
was detected with an outflow mass of about 15 M Q . 
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1 INTRODUCTION 

Massive stars play an important role in the evolution of the interstellar medium (ISM) and galaxies; nev- 
ertheless their formation process is still poorly understood because of large distances, high extinction, and 
short timescales of critical evolutionary phases. In addition, massive stars do not form in isolation but often 
in clusters and associations, which make the environment of massive star formation regions more complex. 

The 6.7 GHz transition of methanol has been found to be a particularly useful signpost to trace massive 
star formation (Minier et al. 2003 Xu et al. 120031 . On the other hand, the maser phase encompasses the 
outflow phase (Xu et al. 2006), which give us another powerful tool to study the dynamics of massive star 
formations. 

IRAS 22506+5944, with an infrared luminosity of 1.5xl0 4 L o , belongs to the Cepheus molecular cloud 
complex. Harju et al. (Q993) made a NH3 map of this region, and found a NH3 core is coincident with the 
peak of the IRAS source. Both H 2 maser (Wouterloot & Walmsley 1986) and SiO (Harju et al. |1998T > 
have been detected. Although searches for 6.7 GHz methanol maser (Szymczak et al. 2000) show negative 
results, recently, we found a weak 6.7 GHz methanol maser in this region. Despite its high luminosity 
and FIR color characteristics of the ultra-compact HII region, no radio emission was detected (Molinari et 
al. [1998). The distances used in literatures for this source range from 5.0 kpc to 5.7 kpc. Here we use the 
value of 5.0 kpc. 

In this paper, we present a multi-line study of this star-forming region. In Sect. 2, we describe our 
observations. The results are given in Sect. 3. We give analysis and discussion in Sect.4, and summarize in 
Sect.5. 

2 OBSERVATIONS 

2.1 The Effelsberg 100 m Telescope 

Observations of the methanol (CH3OH) maser were made using the Effelsberg 100 m telescope in February 
2006. The rest frequency adopted for the 5 1-60 A + transition was 6668.519 MHz (Breckenridge & Kukolich 
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Table 1 Observation Parameters 



translation 


Vrest 


HPBW 


Bandwidth 


Av,-„ 


lcr rsm" 




(GHz) 


(") 


(MHz) 


(km s" 1 ) 


(K) 


i2 C0 J = 1-0 


115.271204 


58 


145 


0.37 


0.10 


13 CO J= 1-0 


110.201353 


61 


43 


0.11 


0.10 


C 18 J= 1-0 


109.782182 


62 


43 


0.12 


0.09 


HCO + J = 1-0 


89.188521 


75 


43 


0.16 


0.10 



" typical value in the scale of T* R . 

1995). The spectrometer was configured to have a 10 MHz bandwidth with 4096 channels yielding a spectral 
resolution of 0. 1 1 km s _1 and a velocity coverage of 450 km s" 1 . The half-power beam width was ~ 2' and 
the telescope has an rms pointing error of 10". The observations were made in position switched mode. 
The system temperature was typically around 35 K during our observations. The flux density scale was 
determined by observations of NGC7027 (Ott et al. |1994| l. The absolute calibration for flux density is 
estimated to be accurate to ~ 10%. The integration time on source was 10 minutes, with a rms noise level of 
~ 0.05 Jy in the spectra. The pointing position was R.A.(J2000) 22 /, 52 m 36.9 ! , DEC.(J2000) = +60°00'48". 

2.2 The PMO 13.7 m Telescope at Delingha 

The 12 CO, 13 CO, C ls O and HCO + J = 1-0 maps were observed with the PMO 13.7 m millimeter-wave 
telescope at Delingha, China, during 2008 November. A cooled SIS receiver was employed, and system 
temperatures was ~ 250 K during the observations. Three AOS (acousto-optical spectrometer) were used to 
measure the J = 1-0 transitions of 12 CO, 13 CO, C I8 and the FFTS (Fast Fourier Transform Spectrometer) 
were used to measure the HCO + J = 1-0 lines. All the observations were performed in position switch mode. 
The pointing and tracking accuracy was better than 10". The obtained spectra were calibrated in the scale 
of antenna temperature T^ during the observation, corrected for atmospheric and ohmic loss by the standard 
chopper wheel method. The grid spacings of the mapping observations were 30". Tabl^T] summarizes the 
basic information about our observations, including: the transitions, the center rest frequencies v rest , the 
half-power beam widths (HPBWs), the bandwidths, the equivalent velocity resolutions (Av res ), and the 
typical rms levels of measured spectra. All of the spectral data were transformed from the T^ to the main 
beam brightness temperature T* MB scale. The absolute calibration for intensity was about 10%. 
The GILDAS software package (CLASS & GREG) was used for the data reduction. 

3 RESULTS AND DISCUSSION 
3.1 spectra 

3.1.1 6.7 GHz CH^OH maser spectrum 

The spectrum of the CH3OH maser detected in this region is shown in Fig. Q] There are two features that 
are separated by about 2.2 km s _1 . The stronger feature is at the LSR (local standard of rest) velocity of 
-53.7 km s , with a flux density of 0.52 Jy, while the other is only about 0.2 Jy. In order to get high signal 
to noise spectra, we did not attempt to refine the position and just integrated the time at the same position. 
Hence, the actual position could be off by 1 arcminute. 

3.1.2 n CO, ]3 CO, C 18 <9 and HCO + spectra 

Spectra of 12 CO, 13 CO, C 18 and HCO + are presented in Fig. [2] The spectra in left panel come from 
the peak of core A (dominant core in Fig. [3]). Both 12 CO and HCO + show remarkable broad line wings, 
with a FW (full width) of 24 km s _1 and 6 km s _1 at lcr level, respectively. Spectra in the right panel are 
correspondent to the peak of core B, and spectra at the conjunctive point of the two cores are given in middle 
panel. Details of the line in positions of the peak, including the line central velocities, the fitted line widths, 
the bright temperatures and integrated intensity were listed in Table|2] 
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Fig. 1 Spectrum of the 6.7-GHz CH3OH maser. The spectral resolution is 0.1 1 km s 



Table 2 Result of molecular line measurements. 



Translation 
(GHz) 


Vlsr 
(km s* 1 ) 


Av rel 
(kms" 1 ) 


T* 

1 MB 
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J T « dv 
(Kkms -1 ) 


'-CO J = 1-0" 


-51.4 


4.3 


22.3 


98.3 


12 CO J= 1-0'' 


-51.5 


3.3 


14.1 


48.9 


13 CO J= 1-0° 


-51.4 


2.3 


9.5 


22.1 


13 CO J= 1-0* 


-51.5 


1.8 


6.1 


11.4 


C 18 J= 1-0° 


-51.6 


1.8 


1.2 


2.3 


C 18 J= 1-0* 


-51.5 


1.1 


0.9 


1.0 


HCO+ J= 1-0" 


-51.1 


3.8 


2.1 


6.9 


HCO + J= 1-0* 


-51.6 


2.0 


0.4 


0.9 



a and b indicate Core A and Core B. 



3.2 Mapping 



3.2.1 n CO, C n OandHCO + 



maps 



Contour maps of the total integrated 13 CO J = 1-0, C ls O J = 1-0 and HCO + J = 1-0 line emissions were 
presented in Figure[3] We used MSX Eband (21yum) image as background images of the integrated contours 
to compare the distributions between gas and dust. The filled triangle denotes 3 millimeter continuum 
peak (Su et al. l2004t . H 2 (Wooterloot & Walmsley [19861 . SiO (Harju et al. [19981 and CH 3 OH masers 
were indicated with the open triangle, star and square, respectively. The IRAS error ellipse is also marked. 
Contour levels are 20% to 90% by steps of 10% of the peak emission with the exception of HCO + J = 1-0 
line, whose contour levels are 10%, 15% and 20% to 90% by steps of 10%. 

From Figure[3] we see that both molecular line and dust emission peak are roughly coincident with the 
IRAS source. 13 CO and HCO + are dominated with a single core. 13 CO shows a little elongation. C I8 map 
clearly shows two cores (Core A and Core B), indicating that such optical thin line traces the inner part 
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Fig. 2 Left panel: spectra at the C O east peak. Middle panel: spectra at the conjunctive point 
of the two C 18 cores. Right panel: spectra at the C ls O west peak. The horizontal dot line is lcr 
level of each line (TableQ}. 



of a molecular cloud than other two lines. The center of the Core A with an angular extent of (70", 60"), 
coincides with IRAS 22506+5944 and MSX peak, indicates that they may be the same source. The Core B 
has an offset of (100", 60") at the north-west of the Core A. The size of Core A is slightly larger than the 
Core B (70", 60"). In order to show the kinematic relation of the two cores we also give the channel maps 
of the 13 CO J = 1-0 lines in Figure H 

We derive the physical parameters of the cores, assuming LTE (local thermodynamic equilibrium) and 
with an abundance ratio [H2]/[ l2 CO] = 10 4 . Given a distance of the source to galactic center, Dec ~ H-4 
kpc ,we adopt an abundance ratio [ 12 CO]/[C n O] - 707 and [ n CO]/["CO] - 93 estimated from the 
relationship [ 16 ]/[ 18 Q ] = (58.8 ± 1 1.8)D GC + (37.1 + 82.6) and [ 12 C]/[ 13 C] = (7.5 + 1.9)D GC + (7.6 + 12.9) 
(Wilson & Rood 119941 1. Excitation temperature is calculated Using equationQ] assuming 12 CO J = 1-0 lines 
are optical thick: 



T* = 5.532 -hn 



5.532 



(r*( 12 CO) + 0.819) 

I3 CO and C 18 J = 1-0 line optical depths, r, are estimated with formulas below: 



(1) 
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Table 3 Physical parameters of Core A and Core B 



Name 


Aa 


A6 


/ 


Tex 


Av" 


N( 1J CO) 
(cm" 2 ) 


N(C ls O) 
(cm" 2 ) 


N(H 2 ) o 


M( Lr£ ) 




(arcsec) 


(arcsec) 


(pc) 


(K) 


km s~' 


(cm- 2 ) 


(M G ) 


core A 


90 


70 


0.9 


26 


1.5 


2.7E+16 


1.7E+15 


2.7E+22 


228 


core B 


85 


60 


0.5 


18 


1.0 


10.7E+15 


7.5E+14 


1.1E+22 


35 



a: Av has been corrected using ^J"" = y Mz/j^^jlf jj] considering line broadening due to optical depth. 
b: H2 column densities were derived using 13 CO column densities, assuming [H2]/[ l3 CO] = 9.3x10 s . 



. , T*{ u CO) 
r( n CO) ~ - ln[l - -5— ] (2) 

T* R o 2 coy 



t(C 18 0) * -ln[l - *;,. ] (3) 



r R (c l *o) 

r R e 2 co) 



13 CO and C ls O column densities are derived using equation|4](Kawamura et al. 1 19981 1 and equation[5](Sato 
et al.[1994). t and Ay are optical depth and intrinsic line width: 

"( C °y- 2 42X 10 1 - exp (-5.29/7-,,) ■ W 

v 1 - exp(-5.27/7V v ) 

The nominal core size, Z, is determined by de-convolving the telescope beam, using equation|6] 

l = D(A 1/2 -(P MB ) m , (6) 

where D is the distance (5.0 kpc), A 1/2 is the area within the contour at the half- integrated intensity of 
the peak and Omb is the main beam size (see TableQ]). 

Core masses are computed with equation|7] where m is the mass of the hydrogen molecule, /j. the ratio 
of total gas mass to hydrogen mass, /j * 1.36 (Hildebrand 119831 1, Nu 2 the column density of H2 and I the 
de-convolved half power size defined above. 

M LTE = M mN Hl l 2 /4 (7) 
The physical parameters derived are tabulated in Table [3] 

3.2.2 Outflows 

Molecular outflows are an important signature of the earlier stage in star formation. An outflow has been 
detected (Wu et al. 2005 ) using the 12 CO J = 2-1 line. A comparison of two different transitions will be 
helpful to our better understanding of the physical properties of outflows. In Fig. [3] we present a similar 
work with the I2 CO (1-0) line. The red and blue lobes are largely overlapped, while the IRAS source is 
located at the center of the outflow, probably the driving source of the outflow. Morphology of the 12 CO J 
= 1-0 outflow is similar to that of the 12 CO J = 2-1 outflow (Wu et al. 120051 , but the former extends a larger 
area than the latter, spreading from Core A to Core B. 

The outflow parameters, except for 12 CO column density which is derived from Snell et al. ( 1988), 
are estimated with the method of Beuther et al. ( 2002). We assume that the gas is in LTE and the line 
wings are optically-thin. Excitation temperature and [H2]/[ l2 CO] abundance ratio adopted are the same as 
Sect. 13.2.11 In order to better define the kinematics of the high gas, we divided the wings into low velocity 
and high velocity segments. The physical properties, including velocity range, size, column density, mass, 
momentum and kinetic energy are summarized in Table [4] 
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Table 4 Outflow properties 



Compoent 


v 

(km s" 1 ) 


Size" 
(pc) 


N(H 2 ) 
(cm 2 ) 


Mass 
M 


P 

(M kms-') 


E* 
(erg) 


red lobe(L) 


(.48.2 -44.0) 


0.8 


1.3E+20 


5.3 


72 


9.6E+45 


red lobe(H) 


(-44.0 -38.0) 


0.8 


4.2E+19 


2.1 


29 


3.8E+45 


blue lobe(L) 


(-54.7 -59.0) 


1.0 


1.0E+20 


6.6 


68 


7.1E+45 


blue lobe(H) 


(-59.0 -62.0) 


1.0 


1.7E+19 


1.0 


11 


1.1E+45 


total 








15.0 


180 


2.2E+46 



a: size of lobes are computed using formula[6] 

Following the method of Beuther et al. (2002), we obtain the characteristic time scale, t m 8.1x 10 4 yr, 
the mass loss rate, M out w 1.8 x 1O~ 4 M yr~' , the mechanical force, F m « 2.2 x 10~ 3 M km s _1 yr -1 , and 
the mechanical luminosity, L„, « 2.2 L . The mass and kinetic energy of the outflow are significantly larger 
than typical values from low-mass star forming regions (Bontemps et al. 1996). 

3.3 evolutionary scenario 

HCO + usually traces the geometrically thick envelope of a core, while C ls O is expected to trace the inner 
part of the core. The C 18 map clearly shows two cores, Core A and Core B, which likely consist of two 
star forming regions. Core A has obvious star forming evidences, such as strong middle and far infrared 
emission, masers and outflows, while Core B is only associated with some cold molecular lines, indicating 
that core A and core B are in different evolutionary stages, Core A at the phase of protostar core, while Core 
B probably at the phase of pre-stellar core. The mass of Core A is more than 200 M , while the IRAS source 
has a luminosity of 1.5 xlO 4 L . According to the relation between mass and luminosity, L ~ M , the core 
mass is around one order of magnitude larger than that of the IRAS source. This indicates there are other 
sources within the core, which are not detected due to the resolution limit of the employed telescope. With 
a rising steep spectrum, the IRAS colors imply that the source is deeply embedded in a dense molecular 
cloud. This IRAS source could be the exciting source of the H2O, SiO and CH3OH masers. Core B has a 
mass of about 35 M , which might form several low or/and medium mass stars in the future. The energetic 
outflow driven by IRAS 22506+5944 covers the whole region, including both Core A and Core B, and 
could greatly affect its surrounding and accelerate Core B to form stars. In summary, the whole region is a 
star formation complex, in which stars at different evolutionary stages live in the same cluster and interact 
with each other. 

4 SUMMARY 

Our multi-line study reveals a star formation complex around IRAS 22506+5944, in which a weak 6.7 GHz 
CH3OH maser was detected. Multi-line Maps reveal a two-core structure: core A with mass of ~230 M 
contains the IRAS source which is driving an energetic bipolar outflow, while Core B, significantly smaller 
than core A, has a large offset from the IRAS source. The two cores are at different evolutionary stages. The 
energy released by the more evolved core (Core A) is influencing the relatively less evolved core (Core B) 
to accelerate its step to form stars. 
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Fig. 3 upper left: contour map of the total integrated 13 CO J = 1-0 line emission in the velocity 
range from -53.4 to -49.5 km s _1 overlaid on MSX E band 21 /jm image, upper right: contour 
map of the total integrated C is O J = 1-0 line emission in the velocity range from -52.1 to -50.8 
kms" 1 . lowerleft: contour map of the total integrated HCO + J = 1-0 line emission in the velocity 
range from -52.6 to -49.5 km s _1 . lower right: Contour map for the 12 CO J = 1-0 outflow. The 
blue wing (solid line) emission was integrated over -60 to -54 km s _I and -49 to -42 km s _1 for 
red wing (dashed line),respectively. Contour levels of the plots are all from 20% to 90% by steps 
of 10% of each peak emission with the exception of HCO + J = 1-0 lines, whose contour levels 
are 10%, 15% and 20% to 90% by steps of 10%. 50% contour levels used to determine core size 
are plotted with thicker lines. The small crosses in the contour plots show the measured positions 
and the ellipses mark IRAS error ellipse. The filled triangles denote 3 mm continuum peak (Su 
et al 2004). H 2 (Wooterloot & Walmsley 1986), SiO (Harju et al 1998) and CH 3 OH masers are 
indicated with open symbols of triangle, star and square, respectively. 
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Fig. 4 Channel maps of 13 CO J = 1-0 lines with contour levels starting at 0.4 K km s _1 and 
separated by 0.3 K km s _1 . 



